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We use the Kubo-Landauer formalism to compute the longitudinal (two-terminal) conductance 
of a two dimensional electron system placed in a strong perpendicular magnetic field, and subjected 
to periodic modulations and/or disorder potentials. The scattering problem is recast as a set of 
inhomogeneous, coupled linear equations, allowing us to find the transmission probabilities from a 
finite-size system computation; the results are exact for non-interacting electrons. Our method fully 
accounts for the effects of the disorder and the periodic modulation, irrespective of their relative 
strength, as long as Landau level mixing is negligible. In particular, we focus on the interplay 
between the effects of the periodic modulation and those of the disorder. This appears to be the 
relevant regime to understand recent experiments [S. Melinte et al, Phys. Rev. Lett. 92, 036802 
(2004)], and our numerical results are in qualitative agreement with these experimental results. The 
numerical techniques we develop can be generalized straightforwardly to many-terminal geometries, 
as well as other multi-channel scattering problems. 

PACS numbers: 73.43.Cd 
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I. INTRODUCTION 



A significant amount of research has been devoted to 
the study of the effect of a periodic potential on a two 
dimensional electron system (2DES) placed in large per- 
pendicular magnetic fields. On the theory side, the so- 
called Hofstadter butterfly^ - the fractal energy spectrum 
associated with the periodically modulated 2DES, in the 
absence of disorder - had been predicted and studied even 
before Hofstadter's seminal paperiSiii^ Later, the trans- 
port properties of such systems were studied by Sti-eda, 
MacDonald and othersiSiSiS After the discovery of the 
Integer Quantum Hall Effect (IQHE)fiii experimental re- 
sults started to become availableiii*i^ Pfannkuche and 
Gerhardts put forward a detailed theoretical treatment of 
transport properties, including disorder effects Their 
theory is analogous to the self-consistent Born approxi- 
mation (SCBA),^"' and it suggests that the splitting of 
one single Landau band into several subbands by the 
periodic modulation can in principle be observed from 
measurements of the longitudinal conductance axx ■ The 
effect of disorder on the fractal structure was found to be 
similar to its effect on the Landau level (LL) structure 
(responsible for the appearance of IQHE): subbands in 
each Landau level are broadened by disorder, but energy 
gaps or pseudo gaps are still open if the disorder is small 
compared to the amplitude of the periodic potential. As 
the Fermi energy sweeps through a subband, the longitu- 
dinal conductance has a maximum and the Hall conduc- 
tance shows a staircase- like jump if the Chern numbers 
in the neighboring gaps are different. Thus, the Hall 
conductance is expected to follow a nontrivial sequence 
of integer multiples of e^/hr^ whereas the longitudinal 
conductance has a series of peaks and valleys as the Fermi 



level moves through different subbands and gaps. 

A recent experiments^ on a high-quality, periodically 
modulated 2DES shows interesting new features in the 
longitudinal conductance, although the periodic poten- 
tial is too weak to produce well-separated subbands (or, 
equivalently, disorder is strong enough to fill in all sub- 
gaps in the fractal structure of each Landau level). Even 
in this case, the longitudinal conductance exhibits re- 
producible oscillatory features in the presence of the 
weak periodic modulation, instead of the single smooth 
Lorentz peak of the unpatterned samples. To our knowl- 
edge, this regime of strong disorder and weak periodic 
modulation has not been investigated in the literature, 
and therefore these recent experimental results do require 
theoretical interpretation. In recent work,^^ we analyzed 
the spectrum and nature (localized or extended) of elec- 
tronic states in such a regime, and showed that simple 
arguments based on these results provide a qualitative 
explanation of the experimental observations. 

In this study, we present a numerical calculation of 
the longitudinal conductance based on models appropri- 
ate for the type of samples used in the experiment of 
Ref. Our model includes a disorder potential and a 
periodic potential with either square or triangular sym- 
metry, with arbitrary relative strengths. In this work we 
assume that these potentials are small enough that Lan- 
dau level mixing is negligible, although the methods we 
develop can be trivially generalized to take such mixing 
into account. Unlike previous theoretical studies dealing 
with disorder effects in QHE, which performed average 
over disorder at the onset of the calculation so that all 
computed response functions are disorder-averaged, we 
calculate the longitudinal conductance from first princi- 
ples for a given disorder realization. This is necessary 
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because the features observed in this experiments^ are 
beheved to be sample-specific. Our calculation is based 
on the Kubo-Landauer formalism^ Our method is in 
principle valid for finite systems with any type of disor- 
der and/or periodic potentials, although computational 
times vary with the sample size and degree of sparseness 
of the Hamiltonian. Regrettably, we have no detailed mi- 
croscopic knowledge of the disorder present in these sam- 
ples and the details regarding their connections to the ex- 
ternal leads; this prevents us from performing meaningful 
quantitative comparisons with the experiment. However, 
the results we obtain for different realizations of disorder 
allow us to qualitatively explain the physics responsible 
for the new features in the longitudinal conductance, and 
to reinforce the arguments offered in our previous work.^^ 
The paper is organized as follows: in Section ^ we 
present the method used for the calculation and the rel- 
evant theoretical considerations. The numerical results 
are presented in Section IIIII while Section IIVI contains 
our conclusions and discussions. 



where I 



is the magnetic length and Hn{x) are 
Hermite polynomials. In the rest of the paper, we con- 
centrate on one of the spin-polarized lowest Landau levels 
(LLL), therefore we set n — 0. This is justified because 
in the experiment both the disorder and the periodic po- 
tentials are estimated to be much smaller than the cy- 
clotron energy hu! and the Zeeman splittingj^^ so that 
Landau level mixing can be safely ignored."'^^'"'^^ Imposing 
cyclic boundary condition in the y-direction leads to the 
restriction Xj = j2TrP/Ly, j = 1, • • • ,N. Xj, the guid- 
ing center, characterizes the location at which individual 
basis states are centered along a:-axis [see Eq. Since 
Xj can vary between and L^, the degeneracy of each 
spin-polarized LLL is iV = L^Ly/ {2ttP). 

Let us define cj to be the creation operator for an elec- 
tron in the LLL: cj|0) = \Xj) (the indexes n = and 
a will be suppressed from now on). In the absence of 
Landau level mixing, the Hamiltonian Tig projected on 
the subspace of the spin-polarized LLL becomes: 



II. THE NUMERICAL METHOD: 
KUBO-LANDAUER FORMALISM 

A. The model 

We consider a two-dimensional Hall sample of rectan- 
gular shape, of size x Ly, with cyclic boundary con- 
ditions in the y-direction and open boundary condition 
in the x-direction, characteristic of a two-terminal ge- 
ometry. Typical sizes we consider are on the order of 
3fim X 3/im. A large magnetic field B of up to lOT is 
applied in the z-direction, perpendicular to the 2DES. 
A rough estimate shows that the degeneracy N of each 
Landau level (LL) is of order 10'', defining the size of the 
matrix to be diagonalized to be 10^ x 10^. For such large- 
size matrices direct diagonalization is time-consuming, 
therefore we look for alternative approaches with a bet- 
ter scaling behavior for large systems. On the other hand, 
the sample size is still small compared to that of the ex- 
perimental sample, which is about 20fnn x 20/im. As a 
result, care must be taken in interpreting the numerical 
results. In principle it is possible to increase the values 
of Lx and Ly, however serious numerical difficulties arise 
when the system size is much larger than the ones we 
consider (these issues are discussed in section Hi C|) . 

The Hamiltonian of the non-interacting electrons con- 
fined in the 2DES is 

Hs^^^Ip+'AY -lgfiB^-B + Vd{x,y) + Vp{x, y). 
2m \ c / 2 

where Vd and Vp are the disorder and the periodic po- 
tentials, respectively. We use the Landau gauge A — 
(0, Bx, 0) throughout this paper, and the complete set of 
eigenf unctions for the n**^ Landau level: 



(r\n,X,<7) = 



Hn 



A/2"n!0F/ 



Xa 



(1) 



huJr 



Af N N 



{X,\Vd+Vp\X,)-c\c 



.(2) 

This looks like a one-dimensional(lD) hopping Hamilto- 
nian, and this is a very appropriate comparison if one 
keeps in mind that the state \Xi) is indeed localized 
within a distance I oi Xi = i2TTp/Ly. 

In order to calculate the matrix elements for the dis- 
order and the periodic potentials, we use the identity 



{X,\e^^-\X, 



^Xi,Xj-qyl 



262 



(3) 



where Q — (q^ + Qy). (The generalization for higher 
Landau levels and/or Landau level mixing is straightfor- 
ward, see for instance Ref. 0|). Let us now consider 
each type of potential separately. 

The periodic potential can be expanded as: 



Vp{r) 



(4) 



where Vg = V*^ because V{r) is real, and {g} are the re- 
ciprocal vectors associated with the Bravais lattice. For 
a square potential, we use Vg = A for all four short- 
est reciprocal vectors g = (±27r/a, 0) , (0, ±27r/a), where 
a is lattice constant, and zero otherwise. Higher order 
components can also be included in the same formalism, 
but result in longer computational time and no quali- 
tative changes. Similarly, for a triangular potential we 
define Vg — —A for all six shortest reciprocal vectors 
g 47r/%/3a(±l,0) ,27r/V3a(±l,±%/3), and zero other- 
wise. The minus sign appears here in order to have the 
minima on the sites of the triangular lattice, as explained 
in Ref. 0|. In both cases the projection of g on the y- 
axis is either or This particular orientation allows 

us to treat these two potentials similarly, since it follows 
that both types of periodic potentials only couple a state 
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\Xj) to itself and to \Xj ± 27rZ^/a) [see Eq. Since 
|Xj ± 27r/^/a) must be in the basis considered, 27r/^/a 
must be an integer multiple of iiiP /Ly, i.e. Ly is an in- 
teger multiple of a. This is consistent with the periodic 
boundary conditions along the y-axis. 

We introduce the integer N^. = Ly/a. From the previ- 
ous discussion, it follows that the periodic potential cou- 
ples a state Xj only to itself, and to the states Xj±pf^. As 
a result, we can divide the total N states of the LLL into 
Nc subclasses (henceforth called the conduction chan- 
nels) using the unique decomposition j — iNc + n, where 
i > is an integer, and 1 < n < Nc. The periodic po- 
tential couples only states in the same channel n; these 
states are distributed equidistantly across the sample, be- 
tween the a; = and the x = edges, and can carry 
currents across the sample (hence the name "conduction 
channels" ) . For simplicity, we require that each chan- 
nel has the same overall number of states, i.e. the total 
degeneracy N is an integer multiple of the number of 
channels Nc, N — pNc- (This condition can be easily 
relaxed.) This imposes a constraint Lx = p{2TrP)/a on 
the values allowed for L^- For the typical sample sizes 
we consider, the constraints on and Ly require only 
minimal adjustments. For instance, we use a = 39 nm^ 
and therefore a sample with Lx ~ Ly k. 2> /xm will have 
around Nc = 70 channels, with around p = 150 states per 
channel. For different values of the magnetic field (dif- 
ferent I values) the length Lx can be kept fixed within a 
few nm by slightly adjusting the value of p. 

For later convenience, we re-label the creation opera- 
tors for states in the LLL as ^ c - „, where j = iNc + n, 
1 < n < Nc- It follows that the periodic potential pro- 
jected on the spin-polarized LLL takes the simple form 



We 



p p-J- 



i=l i=l 

where for the square potential, we have: 



Nc) \ a 



and for the triangular potential, we have: 
— —2Ae sV » 7 cos 



r- , « H 



n 1 

V3a2 V Nc 2 



(5) 

(6a) 
(6b) 

(7a) 
(7b) 



From Equations ^ and Q, it is apparent that the pa- 
rameter controlling the band-structure (in the absence of 
disorder) is the ratio (/)/0o = -^/(27rZ2), where 0o = hcje 
is the elementary flux, and = BA is the magnetic flux 
through the unit cell of the periodic potential. ^ = or 
a V3/2 for square or triangular potentials respectively. 



In particular, if (/)/(/'o = <ilv^ where q and p are mutually 
prime integers, the original LL splits into q sub-bandsfii^ 
In the absence of disorder, there is no mixing between 
different channels [see Eq. ©], and the longitudinal cur- 
rent is just a sum of the currents carried across the sample 
through the individual conduction channels. However, 
disorder introduces scattering between different channels. 
First principles modeling of the disorder in real samples is 
a very difficult and numerically intensive problem iSSiSiiSS 
As a result, we generate the disorder potential using two 
simple phenomenological models described in detail in 
Ref. hj. One is a simple addition of random gaussians, 
while the second model attempts to estimate the proper 
energy scale from considerations of the Coulomb attrac- 
tion between electrons and their donors. Both models 
generate smooth disorder potentials, i.e. with a length 
scale of more than 100 nm, which is very large com- 
pared with the typical magnetic length I (I ^ % nm when 
B ^ lOT). Such long- wavelength disorder is believed to 
be dominant in high-quality samples, as the one stud- 
ied in Ref. (l^. The standard deviation is estimated to 
be 2 to 3 meV, much smaller than the typical cyclotron 
energy.^^ We Fourier-decompose the disorder potential, 
and use Eq.lPJ to compute its matrix elements. The pe- 
riodic boundary condition in the y-direction implies that 
the allowed Fourier components are qy = 2-Km/ Ly, which 
introduces matrix elements between any pair of states 
\Xi) and \Xi±m), where m is an arbitrary integer. It fol- 
lows that different conduction channels are now coupled 
by disorder. The small q Fourier components of the disor- 
der potential are very important, since they describe the 
long wavelength features of the disorder potential. High 
values of q, on the other hand, describe short wavelength 
features of the disorder, which are not well captured by 
our simple phenomenological models. As a result, we 
use a cutoff value of |m| < 36 for Fourier components of 
the disorder potential kept. This value is large enough 
to allow basically exact reconstruction of the disorder 
potential (see relevant discussion in Ref. 17]) but also 
small enough so that the Hamiltonian matrix is still very 
sparse. With this cutoff and in the absence of LL mixing, 
the disorder potential has the general form: 



Na P 

n.n' — l i.i' — l 



(8) 



where 



is non-vanishing only for states within a 



distance \{i — i')Nc + {n — n')\ < 36 of each other. 

From Equations ^ and © it follows that the 
total Hamiltonian for the sample is: 



Hs = 



E 

n=l 



,i=l 



p-1 

E 

i=l 



it- 



,Cj+l,Tl 



h.c 



p 

n,n' — l iA' — l 



(9) 
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where the overall energy shift {huc — 9I^bB<j)/2 asso- 
ciated with the LLL is absorbed in a redefined chemi- 
cal potential. This Hamiltonian can be efficiently gen- 
erated and stored as a column compressed sparse ma- 
trix. In principle, we can directly compute the eigenval- 
ues and eigenfunctions of this Hamiltonian, and calcu- 
late the corresponding Thouless number, characterizing 
its longitudinal conductance.^^ However, this is numer- 
ically very time-consuming. Instead, we use the Kubo- 
Landauer^^i^^i^^i^'^i^^i^^ formula for the longitudinal con- 
ductance (details in section lTlB|l which requires the com- 
putation of various transmission coefficients through the 
sample. The main idea is to link the longitudinal con- 
ductance to the total probability that an electron injected 
into the sample at a; = arrives at x — Lx, or vice versa. 

In order to compute these transmission coeflScients, it 
is necessary to connect the sample to external metal- 
lic leads which allow us to inject into and extract elec- 
trons from the sample. We model each external lead as a 
collection of independent, semi-infinite ID tight-binding 
chains, as illustrated in Figure ^ In reality, the leads 
have, of course, higher dimensionality than one. One 
way to simulate this would be to add bonds (hopping) 
between the various ID chains. We do not add these ex- 
tra bonds for the following reason: the eigenstates of any 
lead with complex geometry in the transverse direction 
have the general structure E'fc^n = Cnik) + Un- Here, k 
is a quasimomentum associated with the longitudinal di- 
rection, a good quantum number given the translational 
invariance along this direction, and n is some discrete set 
of quantum numbers characterizing the discrete trans- 
verse modes supported by the particular geometry of the 
lead cross-section. In other words, any perfectly metallic 
higher-dimensional lead reduces to a collection of inde- 
pendent ID leads (or channels), whose dispersions can be 
simulated by an appropriate choice of ID tight-binding 
chainSfSS. The question, then, is how many channels are 
in each lead, and how are they connected to the states in 
the sample. Since we have no detailed knowledge regard- 
ing the leads, and since one hopes that the main features 
of the longitudinal conductance will come from the sam- 
ple itself, not the details of the lead modeling, we choose 
the following very simple solution: we assume that both 
the left and right-side leads have precisely Nc channels, 
and each one of these channels couples identically to one 
of the conduction channels inside the sample. This is the 
simplest model that satisfies several criteria: (i) the leads 
are perfect conductors; (ii) the conductance of the leads is 
not less than the maximum conductance of the sample^ 
(iii) each conduction channel in the sample has equal cou- 
pling to the leads; (iv) transmission and refiection coef- 
ficients can be easily defined and computed. However, 
any other more complex model for the connection of the 
sample to the external leads can be investigated with the 
formalism we develop here. 

Let us index with n = 1, 2, A^c the Nc independent 
ID channels of each semi-infinite lead, and let c| „ be the 
creation operator for an electron at site i of the n**^ ID 



Left leads 
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Sample 
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FIG. 1: A sketch of the model geometry of the Hall sam- 
ple with periodic boundary conditions in the j/-direction, and 
its contact to leads on both x-axis ends. The lower diagram 
shows a detailed view of the left edge. The first A'^c eigenstates 
\X\), 1X2)..., \XNa) near the edge belong to different conduc- 
tion channels. We assume that each such conduction channel 
is attached to external leads at both edges of the sample. 



chain. We distinguish between the left and right leads 
by requiring that i < for the left channels, respectively 
i > p for the right channels. This notation also avoids 
any confusion between these operators, and the operators 
creating electrons in LLL states inside the sample, which 
correspond to 1 < i < p. The spin-index is suppressed 
everywhere, since in the absence of magnetic impurities 
electrons with different spins travel independently. 

The Hamiltonian describing the external leads and 
their coupling to the sample is then: 




oc \ 

[-^(4nC^+l,n +Ci+l,„c[„) 



h.c. 



Here, t and cq are adjustable parameters, whose selec- 
tion is discussed in section III CI Our model of the leads 
is thus similar to those employed for the study of molec- 
ular conductance)22i2i^ where tight-binding chains are 
used to simulate the gold contacts. In fact, we treat our 
Hall sample as a big molecule wired to contacts. As we 
demonstrate in the following, we only need to include 
explicitly 5 sites for each lead on each side, because we 
apply outgoing or incoming boundary conditions chosen 
so as to give the same result as explicit inclusion of semi- 
infinite leads into our calculation. 

The total Hamiltonian for the sample and leads reads: 



n 



00 

E E 

n— 1 i— — oo 
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+ J2 ^».«;i',«'4,«Ci',„' (10) 

— l n,n' — l 

where ti^n — and e^.n — eo for all z < or i > p (i.e., 
along the semi-infinite leads) whereas inside the sample 
these parameters are given by Equations lO and 

B. Formula for longitudinal conductance 

We calculate the longitudinal conductance using the 
formalism derived in Ref. ^28.] , where it is shown that as 
a function of Fermi energy Ep- 

CTo^AEf) = ^^«a^^b|Gf/(z,z',Sf/?i)|2 

ab 

= jTr{th), (lib) 

Here, z {z') are positions in the asymptotic regions of 
the left (right) leads, which are perfect conductors!^ 
Also, Va(b) s-re group velocities in channels a (b) of the 
corresponding left (right) lead. The Green's function 
G^^'^{z, z', Ep/fi) for a continuous system isi^Si 

riR.A ( / Ef\ 

J J dv^dv\ul{v'^)ub{v^)G''-^ ('■''"''T") ^^^^ 

The above equation is simply a projection of the full 
Green's function G^'^{r,r',Ep/h) onto the transverse 
modes Ua,b of the leads. In the asymptotic regions of the 
leads, the conductance becomes independent of z and z' 
in Eq. (|lla|l because there is no dissipation in the leadsi^ 
In our model each lead is composed of independent ID 
chains, each representing a transverse mode (or channel). 
As a result, the subscripts a and b in Eq. (|lla|l should 
be replaced by indexes of different ID chains on the left 
and right sides, respectively. Eq. Ijllb|) is the Landauer 
formula for perfect-conductor leads, ■^'^'^^ which we now 
identify with Eq. I|lla|l by analyzing the structure of the 
Green's function. 

In Ref. 28] , Eq. (|lla|) is derived for a continuous model; 
the derivation can be replicated with minor modifications 
for a discretized system such as ours. The Green's func- 
tion is now labeled G^{i, n; j, m; E/h), where i and j re- 
place the locations z, z' of the continuous model, whereas 
n, m replace the channel indexes a, b. Consider now the 
scattering state 10^^) of an electron with momentum 
k > and total energy E{k) = eo — 2tcos{ka) injected 
along the n*^ channel of the left lead. The scattering 
state is a solution of the Lippman-Schwinger equation: 

Kij) = ICfc) + ^^.^ (^ - -HoMIe)- (13) 



Here, \4>f^ j.) is the incident wave along the left chain, 

defined by (j, m|0^ j.) = <5n,me'*'^ , which is an eigenstate 
with energy E of the Hamiltonian TLq defined by 

Na OO 

n— 1 i— — oo 

In other words, TYq describes Nc non-interacting, infinite 
ID chains which go through the sample. The scattering 
Hamiltonian Ti. — Tio is then similar to the sample Hamil- 
tonian Hs [see Eq. ©J, except that ei^„ ei^„ — eo and 
ti,n ~^ ti,n + t in order to account for the terms included 
in Ti.Q. As a result, the scattering potential Ti. — Tia has 
non-zero matrix elements only inside the sample. 

The separation of the total Hamiltonian in this form is 
very convenient, because the Green's function Gq {E) — 
[E — TYo + «0]^^ and corresponding group velocity Vk for 
the tight-binding model are well-known: 

{z,n\G^{EUm) = <5„.„, '^P^'^'^' ~ (14a) 

IVk 

Vk = y/U^-iE-eo^, (14b) 
E = eo-2fcos(fe). (14c) 

The quasimomentum k along the leads is measured in 
units of l/d, d being the lattice constant of the ID chains. 
(Since we assume that all ID chains are identical, the 
group velocity and dispersion are the same for all of them. 
Generalization to non-identical channels is straightfor- 
ward, but requires inputs for more parameters). 

One can now prove (see Appendix) that if i < and 
j > p, then G^- (i,n;j,m; I) = (t)+ j^{j,m)e-''''/{ivk). 
On any of the right-side channels, we have (f>n eU' ~ 
tnm{E)e'^^^ for all j > p, where tnm{E) is the ampli- 
tude of probability of transmission from the left chan- 
nel n into the right channel m of an electron with en- 
ergy E. Plugging this asymptotic form of G^ into 
Eq. Hlla|) leads to the Landauer formula, Eq. Ijllb|) . i.e. 

^.AEf) = eyhj:r.,rn\tnmiEF)\'. 

The question, then, is how to efficiently calculate the 
transmission coefficients. In principle, one needs to solve 
the Lippman-Schwinger equation H13|l which, in turn, im- 
plies finding the retarded Green's function G^{E) from 
the Dyson equation (see Appendix). The essential diffi- 
culty for the numerical calculation is to contain the infi- 
nite system in a finite-size scheme of computation with 
appropriate boundary conditions. Our solution to this 
problem allows us to find the scattering solution in a 
very elegant and economic way, which also avoids dif- 
ficulties related to choosing the "proper" value for the 
small imaginary parameter in the denominator. 

To illustrate our solution, assume first, for simplicity, 
that there is only one channel in both the left, and the 
right-side leads. Consider the solution of the following 
inhomogeneous system of linear equations: 
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where sites —M + 1, — Af + 2, • • • ,0 are on the left-side 
lead and sites 1, • • • ,p+M are on the right-side lead, 
and E = cq — 2icosfc, with fc > 0. The matrix in this 
equation is simply the matrix {i\E — H\j), dressed by the 
first and last lines and columns. Along the lead sites the 
matrix has the simple form of a hopping Hamiltonian, 
whereas inside the sample, i.e. for sites between 1 and 
p, it contains the full Hamiltonian of the sample. One 
can now check, using iterations and starting from site 
—M + 1, that the general solution on the left-side lead 
is of the form 0(n) = e**^" + (t){-M + i)e-'HM-i+n) 
all — Af + 1 < n < 0. In other words, it consists of 
an incoming wave with unit amplitude, and an outgoing 
wave with amplitude |(/)(— Af -I- 1)|. On the other hand, 
starting from the site p + M, one finds that the general 
solution on the right-side lead is of out-going type (j){n) — 
4>{M -I- p)e'*=("-*^-P) for aU n > p. Since the matrix is 
non-singular, this inhomogeneous system has a unique 
solution which fixes the values of (j>{—M + 1) and (j){M + 
p). As a result, we can identify directly the reflection 
coefficient for the given energy E as r{E) = (t){—M + 1) 
and the transmission coefficient t{E) = (f>{p + M). 

This approach has a number of considerable advan- 
tages. First of all, the number of sites M kept on the 
left and right-side leads is inconsequential, as long as 
it is greater than 2; the solution obtained through the 
matching of the wavefunction inside the sample is pre- 
cisely the same as for semi-infinite leads. In our simu- 
lations, we keep 5 sites for each ID channel (Af — 5). 
The equivalence of the solution of this finite system with 
that of a system with semi-infinite leads also implies that 
while we work numerically with finite matrices, we have 
a continuous energy spectrum E{k) = eg ~ 2tcosk for 
the electrons injected into the leads (generally, finite sys- 
tems have a discrete spectrum) . Thirdly, there is no need 
to introduce an infinitesimally small parameter iO as in 
Eq. (|13|l . since the first and last line of the matrix in- 
sures that it is no longer singular. Choosing the right 
value for the small imaginary part is always difficult in 
numerical calculations and has to be done very carefully 
in the case of a transfer matrix. ■^^ With our approach, 
we take the limit 77 ^ trivially and obtain the exact 
numerical solution. Finally, one can use various entries 
on the right-hand side of the inhomogeneous system. In 
our calculations, we actually use a more efhcient formu- 
lation, in which all coefficients on the right-hand side. 



except the third entry, are zero, i.e. the equation is 
{E-n)(l) = B, where = (0, 0, 1, 0, 0, • • • .). One 
can verify that this choice injects on the left lead, to 
the right of the 3'^'^ site, an incoming wave with ampli- 
tude l/(2ii sinfc) plus an outgoing wave, meaning that in 
this case, t{E) = 2itsmk ■ <j){M + p). In this case (see 
Appendix) (/)(n) = G^{~M + l,n,E), i.e. some of the 
matrix elements of the retarded Green's function are pro- 
duced. Given the identity of Equations Hlla|) and ljllb|) . 
the two formulations are equivalent, but the second is 
numerically more efficient. 

The generalization to the case with Nc left and Nc 
right channels is straightforward. As already stated, we 
keep only 5 sites for each of the left and right-side ID 
channels, implying that —4 < i < p + 5. As a result the 
total dimension of the matrix A = E — H is Na = pNc + 
5 X 2 X iVc, where NcP = TV is the total number of states 
in the LLL. The matrix elements of A equal the values 
of {i,n\E — H\j,m) for all i,n and j,m values, except 
for (-4, n\A\ - 4, n) = (p + 5, n\A\p + 5,n) = exp {-ik) 
and (-4,7i|^| - 3,7i) = (p + 4, n|A|p + 5, n) = -1 for 
all n = I,-- - , Nc- This insures that proper out-going 
solutions are selected for each channel. 

We then solve uq = 1, - ■ ■ ,Nc systems of inhomo- 
geneous equations of the type J2j,m^'-ri,jm,Xn„ij,m) = 
B„o(i,n), where Bn^{i,n) = 5{n — no)S{i + 2), i. e. cor- 
responds to an electron injected into the rio'^ left-side 
channel. As discussed, the transmission coefficients are 
thentng^m — 2itsmkXng{p + 5,m) (one could choose any 
sites between p + 2 and p + 5), leading to the total lon- 
gitudinal conductance at energy E = eo — 2t cos k to be 

aUE) = j^t^sinHY^ l^no(P + 5,m)|2 (16) 

no — 1 m—1 

Despite having a very large dimension, the matrix A 
is very sparse, and the Nc similar sets of linear equations 
AX — B for a given energy E can be solved very effi- 
ciently using the SuperLU^^ packages. This approach is 
much faster than direct diagonalization and is particu- 
larly well suited for parallelization; as a result, a dense 
grid of energy values E can be investigated. We used a 
cluster of 25 CPUs to scan different energy values in par- 
allel. A typical run lasts for about 10 hours and generates 
5000 data points of axx{E). 
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C. A toy model 

In this section, we analyze the longitudinal conduc- 
tance for a simple toy model. This allows us to under- 
stand the general effect of the lead parameters t and eg 
on longitudinal conductance axx{E). As shown in the 
following section, the shape of the curves axx{E), es- 
pecially for small or no disorder, is rather surprising at 
first sight. It consists of large numbers of very thin reso- 
nant peaks superimposed over a broad maximum. These 
sharp resonance peaks are not numerical errors, and do 
not signal singularities of the type expected in Green's 
functions.'^ (In fact, since we deal with an infinite sys- 
tem, one expects a continuous cut, not individual singu- 
larities, in the Green's functions.) Zooming in the energy 
scale shows that these peaks are features whose width 
scale inversely proportional to the size of the sample (i.e. 
region between the leads). These resonance peaks have 
the same origin as the peaks in the differential conduc- 
tance of molecules attached to metal contacts i.e. 
they correspond to resonant tunneling through the sys- 
tem. However, our 2DES is different from a molecule 
in that it contains about 10^ internal states, hence in a 
small energy interval there can be a large number of such 
resonance peaks. 

A simple toy model that helps in clarifying the ori- 
gin of these resonance peaks is that of an n-bond tight- 
binding model with hopping t2 sandwiched between two 
semi-infinite tight-binding chains with hopping ti (on- 
site energies are all zero for simplicity). Figure [21 shows 
an example of the toy model. 



Left lead 


Sample 


Right lead 


tj 


tl 




tl 


ti 



FIG. 2; A sketch of the toy model. The "sample" with n — 5 
t2 bonds is connected to semi-infinite leads with ti bonds. 

We assume t2 < ii; transmission is then vanishing for 
\E\ > 2t2, since the middle section (the "sample") does 
not support propagating modes at those energies. The 
reflection and transmission coefficients on the two inter- 
faces (dashed line in Fig.j^J can be readily computed, and 
the total transmission rate is calculated either by sum- 
ming up multiply reflected waves, as for a Fabry-Perrot 
interferometer, or by solving the Schrodinger's equation 
directly. The final result is summarized below for a given 
energy E: 

E ~ — 2ti cos(fci) = — 2^2 cos(fc2) 
X = tl sin(fci), y = ^2 sin(A:2), 
a — x + y, P — x — y 

^ ' a4 _^ /?4 - 2a2/32 cos(2nfc2) 




- 



0.2 0.4 0.6 0.8 1 1.2 

E 

FIG. 3: Transmission rate for the toy model corresponding to 
n — 20, tl = 1 and t2 ~ 0.5 and 0.3 respectively. T{E) is a 
symmetric function, so only the interval E > is shown. 



We can see from the last equation that T{E) = 1 
for all energies for which cos(2rifc2) = 1, leading to 
the appearance of a series of thin peaks at the corre- 
sponding energies. The width of each resonance peak is 
roughly 1/n of the band width, i.e. it is inversely propor- 
tional to the system size, as indeed observed in Fig. [S] 
The minima between neighboring peaks correspond to 
T{E) = Ax'^y'^ / [x"^ + V^Y : and are located at energies 
for which cos(2nfc2) = —1. Thus, the larger the system, 
the more and sharper such resonance peaks appear. This 
toy model also provides a check for our numerical scheme, 
both in the single and the multiple-channel cases. 

The decrease in T{E) from its maximum value of 1 
is not due to the resistivity of the "sample" itself, since 
there is no scattering inside the sample in this toy model. 
Instead, it results from the resistance of the contacts, 
which is due to the mismatch ti 7^ t2. (If ti = t2 — > /3 = 
T{E) = 1 at all energies). In fact, if the "sample" 
strip is connected to left and right leads in an infinitely 
smooth manner, the conductance simply becomes e^/Zi 
within the overlap of the spectrum of the "sample" with 
the spectrum of the leads. Similar oscillating behavior 
in conductance due to contact reflection was previously 
studied in Ref. [s^l, where the conductance through a 
short and narrow ballistic channel is calculated exactly. 
The oscillatory behavior is attributed to the "longitudinal 
resonant electron states, the electronic quantum analog of 
the acoustic resonant modes of an open organ pipe" . The 
toy model we introduced above offers a good analogy to 
the continuous ballistic channel. 

In our calculation, we want to concentrate on the 
physics inside the sample, therefore we want to mini- 
mize the additional resistance from the contacts, which 
introduces these extra features in the conductance when 
disorder inside the sample is small. As a result, we have 
to adjust the value of t so that it is close to the magnitude 
of the matrix elements ti^n of Hg. On the other hand, the 
Fermi energy in the calculation is always required to be 
within the spectrum of the leads, Ep G [eq — 2t, eg -I- 2t], 
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so that the leads behave hke perfect conductors. To sat- 
isfy this condition without using a large t value (which 
leads to impedance mismatch), we set the on-site energy 
eo at Ep in each round of calculation, meaning that the 
leads have a "floating" spectrum. The floating spectrum 
is certainly not present in any experiments. We use it 
as a simple way to save computational power. We have 
verified that such "floating" leads do behave as perfect 
metals in the simulation, and that small variations of the 
parameters Cq and t do not change the main features of 
the numerical results. 

The appearance of the large number of resonance peaks 
also illustrates the evolution from small quantum me- 
chanical systems to large macroscopic systems. Quantum 
mechanical quantities like (t{Ef) become rapidly oscillat- 
ing functions as the size of the system is increased. If the 
energy scale for such oscillations is smaller than the reso- 
lution of the measurement or the temperature smearing, 
the measurable physical quantity is a certain average of 
this rapidly changing function within the characteristic 
energy interval. In theory, numerical microscopic calcu- 
lations become inefficient for large systems, because one 
has to sample very many energy values to obtain a good 
picture of the rapidly changing function. In our case, we 
provide curves more suitable for comparison with the ex- 
periment, by convoluting the zero-temperature conduc- 
tance with a sampling function: 



1 



+00 



PT{tJ-,V,e) = — 



1 



1 



1 



(18a) 
.(18b) 



Here, a is the measured conductance at chemical po- 
tential /Lt, temperature T and voltage difference between 
the two edges of the sample V (estimated to be of order 
10~^ eV). This formula is appropriate for non- interacting 
electrons. One can verify that at zero temperature, 



1 



[e{^I + v/2-e)-e{^i-v/2-e)]. 

When y — > 0, pT ^ —dnEp{T, r/)/dr], where nEp{T, r/) 



exp 



/ V-Ef \ 



1 



is the Fermi-Dirac distribution. 



III. NUMERICAL RESULTS 

In this section, we present representative results from 
our calculations. We analyze the interplay between disor- 
der and periodic potentials of various strengths, starting 
with the case of no disorder (pure periodic potential). 

In Fig. 21 the longitudinal conductance of a small sam- 
ple with only square periodic potential is shown. There 
are only 1860 electronic states inside the sample, dis- 
tributed over 30 channels and the magnetic field is such 
that 4'l4>o = llv = 2. Because there is no disorder, differ- 
ent conduction channels do not couple to each other, and 




FIG. 4: Longitudinal conductance of a small sample with 
square periodic potential but no disorder potential. (j)/4>(i = 
q/p = 2, B = 5.44T, = 1.248 /im, Ly = 1.17 nm, corre- 
sponding to a total of 1860 states in the lowest Landau level. 



each of them is similar to the simple toy model discussed 
in the last section. The spectrum of each conduction 
channel (except one) splits into two subbands as expected 
for the Hofstadter butterfly corresponding to g = 2. The 
channel that does not split corresponds to n = 15 in 
Eq. jni, where the cosine vanishes and ei,„ = every- 
where. Each channel contributes to the conductance for 
energies inside its own spectrum, and thus we see many 
sharp resonant peaks on top of a fairly broad conduc- 
tance curve. As already discussed, the sharp peaks are 
due to contact resistance, not the sample itself. The un- 
derlying broad peaks, on the other hand, are a signature 
of the sample behavior. They are simply a reflection of 
the density of states inside the Hofstadter butterfly in the 
clean model. In the vicinity oi E — we see that Uxx is 
bounded by a series of staircases, each of which marks the 
edge of the spectrum of a different conduction channel. 
Thus, the curve for T = (thin black line) is understood 
as a superposition of many ID chains similar to the toy 
model shown in Fig. 13 The two solid curves in Fig.^jare 
the "measured" dxx given by Eq. l(TH)l. corresponding to 
V — 1/ieV and T — 0.1 and 50 mK. At large tempera- 
tures, the resonance peaks are smeared and only the two 
broad peaks are visible, whereas for small temperature 
more detailed features are revealed. 

Figure |31 shows a^x of a sample with only triangular 
periodic potential at q/p — 3. From the semi- log inset, 
we can see gaps between the q — i expected adjacent sub- 
bands. Since the triangular potential is not particle-hole 
symmetric (unlike the square potential) the subbands are 
no longer symmetrically placed with respect to -E = 0. 
In the smaller gap, the conductance is not vanishing, al- 
though it is over 40 orders of magnitude smaller than in- 
side the band. This is a consequence of the fact that the 
sample is finite and tunneling across it is possible even in 
the gap region (although with extremely low probability) . 

FigureElshows the calculated longitudinal conductance 
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of a sample with triangular potential and very small dis- 
order, at q/p = 7/3. In the absence of disorder, we ex- 
pect q = 7 subbands to appear asymmetrically in the 
spectrum. The small disorder closes some of the gaps, 
and only three main subbands are distinguishable in axx ■ 
(Each of the three visible subbands is actually a collection 
of 2 or 3 of the subbands expected in the absence of dis- 
order). The three curves are the same data measured at 
different temperature and different voltage drops across 
the sample. As expected, at low temperature and low 
voltage drop, the measured conductance reveals a vari- 
ety of resonance peaks on top of each broad conductance 
feature. From the semi-log plot in the inset, one can see 
the gap between the left subband and central subband is 
wide open, whereas the smaller gap between the center 
and right subbands is partially filled-in by disorder (the 
value of axx inside the gap is of order 10~^, as opposed 
to 10"''° in the absence of disorder, see Fig. |SJ). This 
proves that narrow subbands are more easily affected and 
therefore more likely to be closed by even small disor- 
der, as one would expect on general grounds. For higher 
temperatures and/or voltage drops, the sharp resonance 
peaks are averaged out and one obtains relatively smooth 
curves with broad peaks reflecting the density of states 
and degree of localization of the sample. 

In contrast to the previous cases. Fig. [7| shows the lon- 
gitudinal conductance for a sample with only disorder 
potential (thick line) and disorder plus a weak triangu- 
lar periodic potential (dashed line). In the disorder-only 
case we see a single broad peak marking the conventional 
integer quantum Hall transition. The curve measured at 
T = ImK and V — IveV is relatively smooth, and its 
flat top at unit conductance indicates the existence of one 
semi-classical orbit extending between the two opposite 
edges. This semi-classical orbit can be seen in Fig. |S1 
where we display the disorder potential used for this cal- 
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FIG. 6: Longitudinal conductance for triangular periodic 
potential plus weak disorder, for q/p = 7/3 (B = 7.33T). 
Lx = 1.998/im, Ly = 2.028/xm. This sample has 7176 states 
divided into 52 channels. The inset shows a semi-log plot of 
the dataset. The amplitude of disorder is about ImeV, which 
is only a fraction of the size of the largest gap. However, some 
of the expected smaller gaps are already filled in by disorder. 





FIG. 5: Longitudinal conductance for pure triangular periodic 
potential with A = 0.01 meV (no disorder) and q/p = 3 
{B = 9.42T). Lx = 2.004/im, Ly = 2.028/im. This sample 



has 9256 states in the LLL, divided into 52 channels, 
inset shows a semi-log plot of the original dataset. 



The 



FIG. 7: Longitudinal conductance for a large disorder po- 
tential. Lx = 2.432fim, Ly = 2A18fj,m, B = 4.71T. 
The thick line corresponds to disorder only, while the thin 
line corresponds to disorder plus a small periodic potential 
{A = O.OlmeV, p/q = 2/3). In both cases T = ImK and 
V — IfieV. The inset is a semi-log plot of the original datasets 
showing the increase of conductance in the off-resonance tun- 
neling regime, induced by the small periodic potential. 



culation. The conduction peak is not centered at E = 
(center of the LLL level) because the disorder potential is 
not fully particle-hole symmetric. However, if one sym- 
metrizes the disorder and averages over many disorder 
realizations, the averaged conductance peak would likely 
be a smooth Lorentzian-type function. The finite (small) 
width is due to the finite size of the sample: all wave- 
functions with a:-axis localization length larger that Lx 
can transport electrons between the leads. 

As expected, the addition of a weak periodic potential 
is not enough to open gaps between the q = 3 subbands 
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FIG. 8: Contour plot of the disorder potential used to calcu- 
late the conductance of Fig. Solid lines are equipotentials 
for meV, and dashed lines for 0.025 meV. These energies 
are located within the central peak of conductance. Parts of 
the contour go along the edge at y — ±Ly/2. 



of the corresponding Hofstadter butterfly; the amplitude 
of the periodic potential is only a small fraction of the 
bandwidth of the disorder-broadened Landau level. How- 
ever, the small periodic potential still has a sizable effect 
on the longitudinal conductance: the magnitude of the 
off-peak conductance has clearly been increased (see the 
inset) especially on the high energy side. Also, the width 
of the central peak is increased by an amount compara- 
ble to the magnitude of the periodic potential A, with 
several more peaks separated by clear valleys appearing 
on both sides of the central peak. 

Figure ini shows another example of the effect of a weak 




FIG. 9: Another example of a weak triangular potential im- 
posed on a large disorder potential. Here, Gxx is plotted as a 
function of filling factor u. Panels (a) and (c) show the results 
with both periodic and disorder potentials, panels (b) and (d) 
are for disorder only. Panels (a) and (b) are T = V = G data, 
panels (c) and (d) are measured at T = 1 mK and V — 1 /iV. 
Parameters are: B = 7.85 T, p/q = 2/5, = 2.432^im, 
Ly = 2.418^m, 11036 states divided in 62 channels. 



FIG. 10: Longitudinal conductance of a sample with a fixed 
disorder potential but varying strength of triangular periodic 
potential. Parameters for this sample are: = 3.107/im, 
Ly = 2.964/im, p/q = 2/5, B = 7.85T, leading to 16872 
states divided amongst 76 channels. The data corresponds to 
T = ImK and V = l^ieV. 



periodic potential on CTj-j.. In this case, we plot a^x 
against the filling factor i'{E) calculated as described in 
Ref. 01 . For disorder only [panel (b) and (d)], we see 
a narrow, double-peaked conductance near half-filling, in 
agreement with the general expectation for the IQHE 
(the double-peak structure is an artifact of the particu- 
lar disorder realization used in this simulation). When 
a small periodic modulation is added [panel (a) and (c)] 
the conductance shows a much more complex shape: the 
central peak is broadened considerably, and several extra 
peaks appear on the low-filling side. 

Finally, Fig. llOlshows the conductance of a disordered 
sample for varying strengths of the periodic potential. 
The disorder potential, which is kept fixed, is plotted in 
Fig. 111! For a strong periodic modulation as compared to 
the disorder {A — 0.25 meV, upper panel), the q — 5 sub- 
bands expected in the Hofstadter butterfly at this mag- 
netic field are beginning to emerge. On the other hand, 
in the disorder-only case (lowest panel) we see a single, 
smooth conductance peak corresponding to the small en- 
ergy interval where percolations through the sample is 
established. Our results show that the interpolation be- 
tween the two cases shows interesting and non-trivial be- 
havior. Although when A — 0.05 meV the disorder is 
large enough to completely erase the butterfly structure, 
the weak periodic modulation has a non-trivial signature 
reflected by large numbers of peaks on both sides of the 
disorder-only main conductance peak. The physical ori- 
gin of these extra peaks has been carefully analyzed in 
Ref. [T^ , where we argued that even a weak periodic po- 
tential can efficiently create supplementary percolation 
paths through the sample for energies within a range 
of the order A from the critical region. Here, we see 
the signatures of these states in the longitudinal conduc- 
tance as well-defined peaks within narrow intervals in en- 



11 




-1.6 -1 -0.5 0.5 1 1.5 



FIG. 11: The disorder potential used in the calculation for 
Fig. 1101 generated by summing random Gaussian scatter- 
ers. Solid and dashed contours are at energies 0.11 meV and 
0.14 meV respectively, marking the extended orbit responsible 
for the conductance peak of the pure disorder case. 



ergy. Even an extremely small periodic potential (such 
as ^ = 0.01 meV) still widens the conductance peak, 
although there are fewer extra peaks visible in this case. 

One technical note: we used a semi-log scale in Fig. 1101 
because the effect of the periodic potential is not so well 
revealed in linear scale. This is partially due to the fact 
that we kept t for the leads constant for all four cases; this 
implies that as A is changed, the impedance mismatch 
increases and the contact resistance becomes more and 
more important, suppressing the value of axx- 



IV. CONCLUSION AND DISCUSSIONS 

In this paper we use the Kubo-Landaucr formalism 
to compute directly the longitudinal conductance of a 
2DES in the presence of periodic modulations as well as 
disorder. The method employed allows us to study indi- 
vidual disorder realizations and thus to analyze sample- 
dependent effects. Our method treats the sample as a 
"big molecule" connected between leads, and is appli- 
cable for any type of one-electron potentials, both dis- 
ordered and periodic. This formulation is particularly 
suitable for long lengthscale, smooth disorder potentials 
such as are believed to be dominant in high quality 2DES, 
because in this case the Hamiltonian is a sparse matrix 
that can be handled numerically very efficiently. Our 
simulations are performed for large mesoscopic samples 
(several microns in linear size) , corresponding to roughly 
lO'' electron states per Landau level. Although this size 
is still smaller than that of most devices used in IQHE 
experiments, this type of calculation can help us under- 
stand the physical processes in detail. This method can 
also be very fruitful for investigating transport in vari- 
ous nano-scale devices, and can also be generalized in a 



straightforward manner for systems connected to more 
than two terminals (leads). 

Here, we concentrate on the interplay between a short 
wavelength periodic potential and long wavelength disor- 
der, and their role in determining the longitudinal con- 
ductance of the 2DES. The phenomenology in the asymp- 
totic limits has been known for a long time. A pure pe- 
riodic potential splits the Landau levels into a number 
of subbands, and the resulting band-structure has frac- 
tal properties as a function of the magnetic field. In this 
case, our simulations show the clearly separated peaks in 
longitudinal conductance corresponding to each subband. 
If a weak disorder potential is added, the smaller gaps 
are closed by disorder, as expected on general grounds. 
Our results are in good qualitative agreement with pre- 
vious studies of the weak-disorder case, using SCBA and 
similar approaches^^*i^ However, unlike our method, the 
SCBA is valid only for small disorder and gives only 
disorder-averaged quantities. On the other hand, our 
method also reproduces the results expected for large 
disorder (no periodic potential) with long lengthscale, 
smooth variation. In this case we find a single, relatively 
smooth peak in the conductance at energies where semi- 
classical percolated orbits connect the opposite edges, 
and negligible conductance elsewhere. 

The new results come when we investigate the cases 
where the periodic modulation is comparable, but smaller 
than the disorder. To our knowledge, this case had not 
been investigated previously. We find that even a weak 
modulation has a non-trivial effect on the conductance, 
with the conduction being significantly enhanced. New 
sharp peaks in conductance develop around the origi- 
nal conductance peak, increasing its width and creat- 
ing complex oscillatory features. Although these peaks 
are clearly due to the periodic modulation, their origin 
is not simply the Hofstadter structure, which is not rel- 
evant for large disorder. The mechanism of en hanced 
conduction was explained in detail by us in Refs. pM7| . 
Basically, the periodic potential helps electrons percolate 
through flat regions in the disorder landscape and thus 
connect localized states to form new conductive states. 
This is in qualitative agreement with recent experimental 
observations^^ which show distinct patterns of peaks and 
valleys in the longitudinal conductance of a periodically 
modulated 2DES, instead of the smooth peak expected 
for unpatterned samples. However, our calculation is lim- 
ited to the band conductance, which has metallic temper- 
ature dependence, and therefore is not suitable to explain 
the temperature-enhanced conduction in the tail of R^x 
seen in the same experiment. We believe that in that re- 
gion, the conduction is due to hopping among localized 
states rather than charge transport through extended 
states. As demonstrated in Ref. ^ small periodic 

potential is very effective in increasing the localization 
lengths and thus the hopping probability. 

Our method of computation has a number of limita- 
tions. One comes from the fact that it is an exact calcu- 
lation at zero temperature for a large system. Especially 
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in the absence of disorder, a^^^ is a rapidly oscillating 
function that requires evaluation at a huge number of 
energy values in order to get a good sampling for tem- 
perature averaging, and this is a major limitation. Sec- 
ondly, the lead modeling is very elementary, although we 
can treat the case of semi-infinite leads exactly. The sim- 
plicity in describing the leads and their coupling to the 
sample comes from our ignorance of their detailed phys- 
ical characteristics. However, one can straightforwardly 
generalize our method to describe more complicated dis- 
persion relations, and different types of channels and/or 
couplings to the sample. In particular, the contact re- 
sistance due to the mismatch between the leads and the 
sample should be minimized as much as possible; this can 
be achieved with a suitable choice of the lead parameters. 

The most obvious limitation of this study is that we 
have little knowledge of the functional form of the disor- 
der and the periodic potential, and how the sample's LL 
wave functions are coupled to the Fermi sea in the leads. 
Until such knowledge is available, detailed quantitative 
comparisons with the experiments are not warranted. We 
use two simple phenomenological disorder models, and 
both give qualitatively similar results. Coupled with our 
understanding of the nature of the wave-function (see 
Ref. this gives us the confidence to claim that such 
effects are genuine and should be observed for smooth, 
long length scale disorder. If the disorder varies on a 
much shorter length scale, as may be the case for quan- 
tum wires and dots (e.g., in Ref. [20ll2^ '). we believe that 
a weak periodic modulation will have a very small, if 
any effect, simply because in this case there are no rel- 
atively flat regions in the disorder landscape where the 
periodic potential plays a dominant role. This is indeed 
confirmed by simulations we performed for models with 
short length scale disorder (not shown here), where ad- 
dition of a weak periodic modulations has no noticeable 
effect on the longitudinal conductance. 

As for the periodic potential, we have been using 
only the components with shortest reciprocal lattice vec- 
tors. However, as experimentalists are designing new de- 
vices with enhanced periodic potentials, the higher or- 
der Fourier components, as well as inter-Landau band 
mixing which has not been considered here, might play 
a role. Such cases can be treated with this formalism, 
the complications being only of numerical nature. As 
more accurate models for the sample disorder and mod- 
ulation, as well as the leads and their coupling, on one 
hand, and more powerful computational facilities, on the 
other hand, become available, this formalism will allow 
for meaningful comparisons with experimental results for 
transport in mesoscopic systems. 
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APPENDIX: SCATTERING STATES 



In the following we demonstrate the relationship be- 
tween G^{—2,no;j,m;^) and 0^„£;(j, w). The scat- 
tering state (j)Xo,E^hT^) = 0';"#,to,B^ Satisfies the 
Lippman-Schwinger equation Eq. (|13|l for n — uq, while 



the retarded Green's function G^{~2,no,j,m,^) satis- 
fies a similar Dyson equation: 



G" = G^ 



1 



E-Ho + iO 



in - Ho)G' 



(A.l) 



where, for simplicity, we do not write the arguments of 
G« and expHcitly Here [E - Hq + iO)-^ = Gl^{E) 
is just the formal expression in matrix form, and iO is 
an infinitesimally small imaginary number which selects 
out-going waves in the analytic calculation. The inhomo- 
geneous term in the Dyson's equation has the form given 
by Eq. (O: 



E 



exp{ik\j + 2\) 



IVk 



(A.2) 



We multiply Eq. (|13|l by and subtract it from 
Dyson's equation Eq. ljA.2(l . to obtain: 

cxp [-ik{2 + j)] - exp [ik{2 + j)] 



+ 



IVk 



E-Ho + iO 



(A.3) 



This 



where ip = G^{~2,no,j,m,^) 
equation can be solved by noticing that Ti. — Ho has 
the same structure as Tis, whereas the inhomogeneous 
term has support only on the no"^ lead, outside the 
sample. It follows that the first order solution = 

^no .m^{ 2 - 

and therefore: 



2sin[^^2+2)l actually the exact solution. 



G^(-2,no,j,m,-) = 



E. _ e^'^^^tfiO''"^) 



IVk 



J) 



2 sin [k{2+j)] 



Vk 



(A.4) 
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This shows that up to a multiplicative constant, the 
Green's function reproduces the scattering state on the 
right-side leads (j > p) , and its asymptotic form there is 

therefore -. — tnn me*'^-' . 
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